Y-junction connecting Luttinger liquids: fixed point structure and conductances 
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We study the transport properties of three Luttinger liquid wires (with possibly different inter- 
action strength), connected through a Y-junction, within the scattering state formalism. We first 
formulate the problem in current algebra language and focus on the case of a symmetric set-up, 
for which the scattering matrix and the matrix of conductances is parametrized by two variables. 
For these we derive coupled RG equations, first in a ladder summation up to infinite order in the 
interaction. The fixed point structure and the implicit solution of these equations is presented. It is 
shown that higher order terms beyond the ladder approximation do not change the scaling behavior 
near the fixed points. For sufficiently strong attractive interaction a new fixed point with unusual 
properties is found. 



I. INTRODUCTION 

Electron transport in strictly one-dimensional quan- 
tum wires is governed by the Coulomb interaction be- 
tween electrons. This is spectacularly demonstrated 
by the fact that within the Tomonaga-Luttinger liquid 
(TLL) model in the limit of temperature T — > the con- 
ductance of a quantum wire with finite barrier tends to 
zero, provided the interaction is repulsive and assumes 
the maximum value G = 1 (in units of the conductance 
quantum Go = e^/h ) for attractive interaction. The 
latter behavior appears independent of the strength of 
the scattering at the barrier, and can be traced back to 
the formation of Friedel oscillations of the charge den- 
sity around the barrier, leading to an infinitely extended 
effective barrier potential in that limit. 

TLL behavior has recently been studied experimen- 
tally in carbon nano-tubes."'^'^ In a future nanoelectron- 
ics constructed out of quantum wires, junctions of three 
or more wires will necessarily be involved, requiring a 
knowledge of the fundamental transport behavior of elec- 
trons in such structures. It is known that such systems 
exhibit rather rich TLL effects which have been the sub- 
ject of a number of recent papers. Much of the work in 
this field has used the bosonization method, which gives 
rise to the problem of how to preserve the fcrmionic char- 
acter of charge carriers. When the number of wires meet- 
ing at a junction exceeds two, the Klein factors, which 
give Fermi statistics to the bosonized operators of differ- 
ent wires, are more difficult to handle.'^ Oshikawa ct al. 
recently introduced^^ a new method to study this prob- 
lem, mapping it into the dissipative Hofstadter model 
(DHM), which describes a single particle moving in a 
uniform magnetic field and a periodic potential in two 
dimensions and coupled to a bath of harmonic oscilla- 
tors. When the three quantum wires enclose a magnetic 
fiux, the mapping to the DHM also allows to identify a 
new low energy chiral fixed point with an asymmetric 



flow of current that is highly sensitive to the sign of the 
flux. The paper by^* is giving a systematic and thor- 
ough treatment of the three lead junction problem, and 
we shall refer to it later, when we compare our results 
with theirs. It should be noted, however, that in spite 
of their success, the above mentioned methods do not in 
general allow to determine the fixed point (FP) structure 
of the theory in an unbiased way, as they require a certain 
knowledge of the cxistcmcc^ of FPs as a starting point. 

In this paper we follow a different route: we de- 
scribe the transport properties in fermionic language, 
thus avoiding the problem of Klein factors and the ne- 
cessity to assume the existence of FPs. As shown in our 
previous work, one may recover the known results on 
transport through a TLL wire with barrier and obtain 
new results not accessible by other methods within the 
fermionic formulation.^^ Moreover, as this method is for- 
mulated in a scattering wave picture, the connection of 
the interacting wire to non-interacting reservoirs is nat- 
urally included. While this formulation has previously 
been regarded as being restricted to weak coupling^'^°'^^ 
we have shown in an earlier work that it may be extended 
to the strong coupling regime in a controlled way by using 
an infinite resummation of perturbation theory (ladder 
summation). The result obtained in this approximation 
is universal, i.e. does not depend on the ultraviolet cutoff 
chosen. 

We restrict ourselves to a symmetric time-reversal in- 
variant setup, which is characterized by two independent 
components of the conductance tensor. This setup in case 
of one wire without interaction is relevant to the problem 
of tunneling into a Luttinger liquid and was studied in^^ 
for weak interaction. There it has been found that the 
asymmetric fixed point A (perfectly conducting wire and 
a vanishing tunneling amplitude into the tip, see below) 
becomes unstable in the case of weak repulsive interac- 
tion once the description of the Y-junction is not reduced 
to the simple tunneling model. A further remarkable re- 
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suit of that work is the finding of a nonmonotonic be- 
havior of the conductance (e.g. as a function of temper- 
ature), in certain cases. We extend that work here to 
arbitrary interaction and find similar behavior at strong 
interaction. 

The resulting coupled RG equations will first be ana- 
lyzed with respect to their fixed point structure. We find 
four fixed points, labelled N,A,M,Q, where at N (Neu- 
mann) the three wires arc totally separated, at A (asym- 
metric) the third wire is separated, while the main wire 
is perfectly conducting, at M (mystery) and Q (quaint) 
the conductances take an intermediate value, depending 
on the interaction strength. All fixed points are located 
on the boundary curve marking the allowed area in the 
plane of the two conductances. For repulsive interaction 
only N is stable, A being a saddle point. For attractive 
interaction N is unstable and M is stable. For stronger at- 
tractive interaction, beyond a threshold value A becomes 
an additional stable fixed point. The conductances are 
calculated analytically for several special cases of inter- 
est. We present and discuss the power law exponents 
appearing in the various regimes. 

In order to check the reliability of the ladder summa- 
tion result, we calculate all scale dependent terms (in the 
limit of T = ) up to third order in the interaction. We 
find several terms additional to the ladder terms. These 
terms are nonuniversal and subleading, in the sense that 
they disappear for the case of repulsive interaction, when 
one of the stable fixed points is approached. For attrac- 
tive interaction, the properties of the new fixed point M 
appear to be affected by the three-loop terms beyond the 
ladder approximation. 

The paper is organized as follows. In section II we 
define the model of a Y- junction considered here. We 
introduce the currents in the SU(3)-representation and 
define the Hamiltonian in the chiral representation. Sec- 
tion III is devoted to an analysis of the scattering matrix 
and the derivation of the conductance tensor in terms 
of the S-matrix components. We consider a simple tun- 
neling case and the totally symmetric case, both charac- 
terized by only a single free parameter in the S-matrix 
and the general Y-junction case with two parameters and 
therefore two independent conductance components. In 
section IV we show how the RG equation for the simple 
tunneling case is extracted from perturbation theory in 
the ladder approximation, and discuss the ensuing con- 
ductance. The RG equation for the general Y-junction 
case is presented in section IV B. We further present the 
results on the nonuniversal terms in third order perturba- 
tion theory in section V. Finally, the RG phase diagram 
is presented and discussed in section VI. 



II. MODEL OF A Y-JUNCTION 

We consider the setup of a Y-junction: a quantum wire 
of finite length 2L connected to noninteracting leads, and 
a tip with attached wire forming a contact (junction) in 



the center of the wire, which for simplicity is assumed to 
have the same length L, also connected to a noninteract- 
ing lead. We choose the origin of the coordinate system 
at the junction and denote the two halves of the main 
wire by indices 1 and 2 and the tip by 3. We consider 
spinlciss fermions. In the following we will refer to the 
so defined symmetric (leads 1, 2) set up as a Y-junction. 
A particularly simple version to be considered below is a 
simple tunneling junction (no tunneling barriers, only on 
site tunneling). Below we will also briefly consider a to- 
tally symmetric junction (with respect to interchanging 
any pair of wires). 

The electrons are assumed to interact via a short- 
ranged interaction of arbitrary strength and sign between 
incoming and outgoing fermions. In the present case, in- 
stead of defining right- and left-movers, it is more ap- 
propriate to speak of incoming and outgoing waves, with 
respect to the junction, or "origin" . 

The short-range interaction is characterized by the am- 
plitude g in the main wire ( 1 and 2) and a different am- 
plitude 53 in the third wire. This interaction takes place 
at X < i, measuring from the origin. In order to spatially 
separate effects of the potential scattering at the origin 
and the interaction, we also assume that the interaction 
takes place at a; > oq, whereas the potential scattering 
happens at a; < oq. The scale ao appears below as short- 
range cutofi^ of the theory. 

The Hamiltonian of the system is 




+2-KVFgji)]^ini)j,in'^],outi'3,out (2) 



where we did not write explicitly the range function of 

the interaction. In the chosen parametrization the am- 
plitudes Qj are dimensionless and enter the subsequent 
expressions in the most convenient form. 

We unfold the setup in the usual way,^'^'^^ by putting 
the incoming fermions on the negative a;-axis and the 
fermions, which have passed through the junction, at 
positive As a result, we have a three-fold multiplet 
of right-going fermions with non-local interaction. This 
procedure is depicted in Fig. 1 

The boundary condition at the origin is described by 
the scattering S-matrix as follows. For elastic scatter- 
ing by the central dot, the outgoing fermions at the ori- 
gin are connected to the incoming ones by the relation 

V'i.out = Slrni>ln,in- We choose S in the symmetric way 
as described at length in the next subsection. 

In the scattered states representation, the right- and 
left- going fermionic densities acquire the form 

i>3,in{x) = tlJj{x), Vj,out(a;) = {S.ll))j{-X) 
V'],i„(a;)V'j,m(a;) = pj{x), Vj>„t(a;)V'j,out(a;) = Pj{-x) 

(3) 

here and below we use the notation A = .A.S. 



3 




(a) 



FIG. 1: (Color online) (a) The geometry of a Y-junction 
is shown together with currents of incoming and outgoing 
fermions, the type of fermions in the scattered states rep- 
resentation is indicated by its color. The local short-range 
interaction between the fermions is shown by a wavy line, 
(b) The equivalent representation in the chiral fermion basis 
after the unfolding procedure; the initially local interaction 
becomes non-local, as described by Eq. (8). 



Let us explain the meaning of these quantities. 
We introduce a muhiplet of incoming fermions ^' — 
("01, 'ip2, V'a)- The incoming density pj = 'i'^ .pj.'i' is given 
by a diagonal matrix, i.e. 



Pi 



P2 = 



P3 




(4) 



'0 

,0 



may now be written. 

Pl{'~x)pi{x) + P2{-x)p2{x) 



(5) 



+ l¥{-x)X3^{-x)-^\x)X3^{x) 



where A_,- — .Xj.S, or simply 



PiPi + P2P2 = 



1 



+ 2'^3A3 = P+P+ + P-P- 



P3P3 



1 



.A, 



1 



fAo 



(6) 



(7) 



3 / VV6 



V3 



V6 



Here we defined p 

and p- = 75-^3- 

In terms of these quantities, the Hamiltonian takes the 
following form (from now on we set = 1) 



/oo 
dx ^l^^V-j 
-00 -1 o 



J = l,-3 (8) 

-I- 27re(a <x < L) {gp+p+ + gp-P- + gsPsPs) 

where p± , ps refer to —x and p± , ps refer to x, and Q{aQ < 
X < L) is equal to 1 within the interval specified and zero 
elsewhere. 

For simplicity we do not consider here the so-called 
(74 part of the fermionic interaction, i.e. the terms 

7rfF5j[('0],„V'j,m)^ + {ipj,out'^j,out)'^]- It is known that 
the (74-interaction can be absorbed into the redefini- 
tion of the group velocity inside the interacting region, 
vpj — w_f(1 For finite L one can show that g4 does 

not lead to a renormalization of d.c. conductance, which 
is the quantity of our interest below. The effect of g^ on 
the a.c. conductance can be analyzed, e.g., following the 
guidelines in^^. 



Here and below Xj, with j = 1,...8 are the traceless 
Gell-Mann matrices, discussed in the Appendix A. In 

addition to these, we use also the matrix Aq = y^lj 
which is proportional to the unit matrix, 1. 

The outgoing densities are given by pj = 
SjkS'^iSkmipjipm = .S^ .pj.S.^, SO that in the 
matrix representation pj = S^.pj.S. We will mostly 
omit the hat sign over p below. 

In terms of the above densities the interaction terms 



III. S-MATRIX AND CONDUCTANCES 

A. S-matrix 

The most general S-matrix is defined as follows 
/ ri, ti2, ti3\ 

S = ^21, ^2, t23 (9) 

\t3i h2 rs J 

where rj is the reflection amplitude for wire j, and tj^ is 
the transmission amplitude between wires j and k. The 
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matrix S is unitary S = 1, which allows its param- 
eterization via the exponential, S = exp X]j=o ^i-^i) • 
Obviously, there is a redundancy in the description by 
the nine real- valued parameters 9j in the exponent, as 
only the densities and not the fermion amplitudes enter 
the observables. 

The number of physically relevant parameters for the 
description of the 5-matrix may be determined as fol- 
lows. One may fix the relative U{1) phase between the 
ingoing and outgoing electrons in each wire, by demand- 
ing that the reflection coefficients are real valued. This 
excludes the A3 and As components. It turns out that it 
is more convenient to keep the Ag component, at the cost 
of introducing some redundancy, see below. In addition 
the overall phase described by the Aq term may be set to 
zero. 

A further reduction in the number of independent pa- 
rameters may be derived from the following considera- 
tion. In the limit of three almost detached wires, we 

have S' ~ 1 -|- «Xlj=o ^i^'^i' ^^'^ ^3 ^- "^^^^ limiting 
case elucidates the meaning of the 6j as the tunneling 
amplitudes between the corresponding wires. For exam- 
ple, the infinitesimal hopping between wires 1 and 3 of 
the form (tiaV'aV'i + h.c.)\x=o leads to 6*4 = Re{ti^)/vp, 
05 = Im{ti3)/vF with vf the Fermi velocity, cf.^*. It 
should be noted that this correspondence between the 
Hamiltonian and the ^-matrix holds only in the limiting 
cases, and starts to depend on the definition of the reg- 
ularization procedure in higher orders of tki, i.e. beyond 
the Born approximation.^^ 

By using the above correspondence of the (small) 0j 
and the tunneling amplitudes it becomes clear that with- 
out loss of generality one may require real- valued hopping 
amplitudes ^13 and ^23, i-e. coinciding phases for hopping 
from the third wire into wires 1 or 2 (in the absence of 
magnetic flux). This makes 6*5 = 6*7 = 0. After that 
the phase of the complex valued amplitude ti2 between 
the wires 1 and 2 is fixed. In the absence of magnetic 
fields it must be zero^^ anyway, which further reduces 
the number of components, since it requires 62 = 0. 

Thus in the time-reversal invariant situation, the S- 
matrix can be parametrized by three angles, Oi,94,9q. 
In the presence of a magnetic field a fourth component, 
62, appears. We do not consider this case here. 

In this paper we concentrate on the case of tunneling 
into the center of a Luttinger liquid wire. This amounts 
to full symmetry between the wires 1 and 2, and it re- 
duces the number of independent parameters further to 
only two, 61 and 64 = 6q. 

As mentioned above, wo will in addition keep the angle 
^8- The reason for the introduction of the Ag component 
is two-fold. First, the explicit analytic expressions for S 
below are somewhat simplified. Second, even if we choose 
to start without Ag component, it will be generated dur- 
ing the renormalization process, as we show below in Sec. 
IV B. This means that the three angles 9j,j = 1, 4, 8, are 
not independent. Below we will identify proper combi- 



nations of these variables forming a minimal set of two 
independent variables, in terms of which all the other 
quantities may be expressed. We then parametrize the 
three angles in terms of a set of new angles 9,(j),ip as 
follows 



V OjXj = 9cos(j){Xi + Xe)/V2 

U (10) 
+ \{esm<i> + i})X^ + ^{39sm(t>-tP)Xs 

Explicitly we have the representation of the S'-matrix 
in terms of three angles 6, (j), tp 



'rx h t2 
S = e'^/^ I h n t2 
J2 t2 r2j 



(11) 



with 



i(e 
i(-e 



tl 

r2 = cos 6 — i sin 9 sin 



cos tl + I sm tl sm q 
^'^ +cos9 + i sin 9 sir 
sm 9 cos 6 



(12) 



To identify one of the variables as redundant, we use 
the transformation {9,4','ij)) — )• (^,V',7), with 



tan 7 ■ 
cos 9 



sin 4> tan 9, 
cos 9 



V' = V' + 7, 



(13) 



cos 7 



After some calculation we find the components of the 
^-matrix (up to an overall phase factor e'^'^'^~'^^/^) as 



n = i(e-''/' cos 61)6^-^ 

= l{^-e-^^^ +cos9)e''< 
t2 = ^sin9 
r2 = cos 9e~^"' 



(14) 



We show in the next subsection that the new third 
variable 7 does not appear in any of the components of 

the conductance. In addition, we will see later in Sec. 
IV B that 7 does not take part in the renormalization 
process. 

It will turn out to be useful to consider the most 
elementary case of tunneling from the third wire into 

the main wire (no next neighbor tunneling or tunnel- 
ing barrier) separately. In that case the S-matrix is 
characterized by a single angle, 9 = \/2^4 = V^9q, 
whereas 9i = 0. In terms of the above variables we have 
{9, </>, ^) = {9, 0, 0) or else {9, V, 7) = {0, 0, 0). 

In the fully symmetric case of a junction of three 
identical wires (identical interaction strength and reflec- 
tion coefficients) we have only one independent angle 
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characterizing the S-matrix, 9 = 9i = 64 = 9^. The 
connection to the above (0, 0, ■0) , {9,ip,j) is not very 
transparent, and up to an overall phase factor we have 
n = r2 = (2 + e3*^)/3, ti=t2 = (-1 + e3'^)/3. 



Conductances 



in the general form 

Gaa Crab 



G 



Gba G, 



bb 



1 - iTr(A3A3 



^Tr(A3AB) 



(18) 



-^Tr(AsA3), I (1 - irr(AsAs))_ 



The observables we concentrate on in this work are 
the linear conductances. We first discuss the question 

of the number of independent linear conductances. The 
matrix of conductances, Gjk, is defined through the re- 
lation connecting the current Ij in a given wire with the 
electric potentials Vk in all the leads as Ij = J2k GjkVu- 
The current in the jth wire is given by Ij{x) = 
evp{{pj^in{'-r)) — {pj^out{—x)))- The electric potentials 
Vk give rise to the following source term in the Hamil- 
tonian, Hy = eJ_^dxJ2k^k{Pk,in{x) + 

Pk,out 

i-x)). 

In linear response theory^^ the currents are given by 

Iji^) = !-r^dy({pj{x)pk{y)) - {Pj(-x)pk{y)))Vk In the 
static limit and in the absence of interaction the response 
functions may be evaluated to give 

Gjk = 5jk - Tr{pjPk) = djk - ISjkl"^ (15) 

One easily verifies that the charge is conserved, Gjk = 
0, and that applying equal potentials to all wires pro- 
duces no current, J2k^jk = 

In view of these conservation laws, it is more instruc- 
tive to discuss the current response to certain bias volt- 
ages. Let us define -^6, -^o) = G.{Va, Vb, Vq), with 

Va^iVi-V2), /„ = (/i - /2)/2, 

Vb = {Vi + V2-2V3)/2 4 = (/i + /2 - 2/3)73, (16) 
Vo = (Vi + V2 + V3)/S, lo = ih +I2 + h)/S, 

In this notation, we seek the retarded response of the 
current Ia{x) = ^evpii'tp^ Xaiiix)) - (ip^XatlJi-x))), to 
the source term , e.g., Hvm = § J*^^ dxVaip^ X3tp{x). We 
can write these combinations symbolically in the static 
limit as 

Hv = 5I4A3 + ^VfeAs -I- ^^0-^0, 

/a = 5(A3-A3), /6 = ^(A8-A8), (17) 

-fo = ^(Ao-Ao)-^O, 

As expected, /o = 0; in addition we may choose the zero 
of the electric potentials such that Vq = 0. It can be 
shown^^ that the conductance in the d.c. limit is pro- 
portional to the trace of the product of current and the 
source vertices. It follows then that the line and the row, 
corresponding to Iq and Vq are identically zero, and we 
omit them for clarity below. The remaining four compo- 
nents are non-zero and we obtain the conductance matrix 



1. Conductance of Y- junction 

For the particular choice of S given in (14) above we 
have 

_ /i[l-cos6icosV'], \ 

We see that the new third variable 7 does not appear 
in any of the components of the conductance. In addi- 
tion, we show below that 7 does not take part in the 
renormalization process. Hence, in line with the above 
argument, the general symmetric case is completely de- 
termined by two independent parameters, which can be 
ultimately chosen as the two non-zero components of the 
conductance matrix, Gaa = Ga^ Gbb = Gfa, Eq. (19). We 
will present two coupled RG equations for the two con- 
ductances Ga,Gb in Sec. IV B below. 

It follows from the parametrization (19) of the con- 
ductances that the physically accessible area in the two- 
dimensional space of conductances is not simply given by 
the unit square, but is defined by 

< G6 < 1 - 4 (G„ - i)' (20) 

We will see below that the boundary of the physical 
regime in GojGfc-space defined by (20) plays a special 
role in that all fixed points of the problem are located 
on the boundary. We will use this fact to our advantage 
when we consider the simplified structure of the RG-fiow 
along the boundary curve. 

In the most general case, the matrix of conductances 
has four independent components, in accordance with the 
analysis of the number of relevant parameters character- 
izing the general S'-matrix given in^^. 

It is worth to note the following property of the bound- 
ary (20). It was proven in^^ that the boundary of the 
region of allowed conductances, defined by Eq. (15), cor- 
responds to the matrices S which can be made entirely 
real by certain "rephasing", i.e. multiplication of rows 
and columns of S by phase factors. It was also proven 
that one can recover S from G up to this "rephasing" .^^ 
In our particular case (14) the "rephasing", leading to 
the real-valued last line and last column, is achieved by 
diag[ie~''^ ,ie~''^ ,l].S.diag[i,i,e^^]. It is then clear, that 
the boundary corresponds either to ijj = 0, vr, or = 0, tt, 
which is indeed the shown below (see also^^). 
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2. Conductance in the simple tunneling case 

Let us now consider the special case when the S-matrix 

is characterized by only a single parameter, 6 = 6^ = 9^. 
This corresponds to the simplest model of tunneling from 
the tip into the wire. Then we have 

A3 = cos 6* A3 - sin^(A4 - Ae) 

As = -^sin^filAi - y|sin26i(A4 + A6) (21) 

l + 3cos26»^ 

+ As 

and the reduced conductance tensor follows as 

We note that the conductances satisfy the relation 
Gft = 1 — 4(Ga ~ ) implying that the simple tun- 
neling case traces the boundary of the physically allowed 
region in Ga, Gfj-space, see (20). As we will see below this 
case is more than a specialization to a very simple phys- 
ical model. It actually already contains the information 
on the fixed point structure of the RG flow of the general 
model. Its advantage is that the fixed points and even 
the conductances may be calculated analytically. The 
stability of the fixed points cannot be decided in this re- 
stricted model, as runaway flow away from the boundary 
may occur (see Sec. IV B below). 

3. Conductance in the fully symmetric case 

For the fully symmetric case with ti = t2 = g(e^*^ — 1) 
we have from Eq. (19) 

which shows that the maximum transparency of the fully 
symmetric junction is reached at 6* = 7r/3 and corre- 
sponds to Ga = 2/3, Gb — 8/9. These values are smaller 
than the maximum (unity) values of the individual con- 
ductances, which according to (20) are attained at the 
points Ga = 1/2, Gb = 1 and Go = 1, = in the 
Ga — Gb - plane. 

IV. RG EQUATION: UNIVERSAL 
CONTRIBUTIONS 

A. Simple tunneling case 

1. Lowest order 

The renormalization of the S'-matrix due to interaction 
can be understood by considering the simplest diagrams 



in perturbation theory. To illustrate our approach, we 
consider first the special case of the S'-matrix, (9, (p, ip) = 
{6, 0, 0) in Eq. (10). Our consideration follows closely our 
analysis of the simpler case of one impurity in a Luttinger 
liquid wire.^^ In comparison with this previous case, we 
now have a few channels of interaction in the Hamiltonian 
(8). For the reader's convenience, we now outline the 
basics of our approach (see also Sec. V). 

Consider the causal Green's function for fermions 
Gij{x,t]y,t') = -i(TtVj,m(a;,i)V'],o«t(-2/.*'))- In the 
non-interacting system {gj = 0) and given our model as- 
sumption of equal Fermi velocities in the wires, we have 
Qij {x, t; y, t') = Q{t' — t,y — x)SJ^, according to (3). Here 
the scalar quantity Q describes the kinetic part of the 
Green's function in the scattering states representation. 
For the Wick-rotated (imaginary time) quantity G we 
have 

g{iuj ; z) = -ivphign{Lu')e{ujz)(r'^'/'''' . 

The renormalization of the junction is obtained by con- 
sidering the d.c. limit of Gimix < —L,y > L,a; — )• -|-0), 
and the contribution of interaction terms order-by-order. 
Without interaction we have in this limit —iS*^, when 
we set vp = I. The renormalized S'-matrix is then de- 
fined by -iSi^lr = limu^+o Gim{i^; X < -L,y > L) The 
first correction to Qim{x,0;y,t) in the basis of scattered 
states is of the form 

SQimix,t;y,t') = 2'Kgk J dz j dT{i()i{x,0) 

X (4(-z,r)pr'>..(-z,T)) 

X (V], {z, r)g'^" V,-. {z, r)) V] {y, t))S*^ , 

where summation over repeated indices is implied, and 
the matrices pk and pk are given in (4). The two possible 
ways of contraction of the fermion operators lead to the 
expressions 

2ngkg{T, -z - x)p'i'g{t' = 0, 2z)p(^g{t -T,y- z)S*^ , 
and 

2Trgkg{T, z - x)gQ{i! = 0, -2z)(}'^'Q{t - T,y + z)S*^ . 

Multiplying these expressions by e"^* and integrating 
over T, we note that the dependence on x and y dis- 
appears in the limit w — >■ 0. The renormalized value of 
—iSi^ is hence of the form 

-isUr = -iSL - 2m / dz [p'i'g{t' = 0,2^)pi'^' 

+ p'i:g{t' = Q,-2z)(^i')s*^, 
or, symbolically, 

S*.S=l- 2mgu I dz (^(0, 2z)pkPk + Q{Q, -2z)pkPk) ■ 
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The above eorrcction may be interpreted as the 
self-energy Fock diagrams Sf, since the renormaUzed 
propagator may be represented generally in terms of 
the self-energy as S as Qr{—L,L) = Qo{—L,L) + 
QQ{—L,y)Yi{y,z)Qo{z,L), where integration over y,z is 
implied. Using G(0,2z) = —i/Anz, we may write this 
correction (first at 53 = 0) as 



dz 

T.F = -g j — {[p+,P+] + [p-,P-]) 
= --A ^ ^ '^^ sing 1 + cos 6') 



(24) 



with A = ln(Z//ao). We see that the matrix Green's func- 
tion receives an off-diagonal static correction, signalling 
the necessity to redefine the rotation angle 9 in the rep- 
resentation (10) of the S-matrix. It is remarkable that 
the above correction may be interpreted as a change of 
the angle B, S6 , as it is directed along the initial vector, 
i.e., along A4 -I- Ag. This means that in the case of simple 
tunneling we are allowed to consider only a single renor- 
malization equation for 9, as opposed to a set of three 
RG equations for 9, tp and 7. 

In the presence of the interaction in the third wire, con- 
tributing a term gsPsPz in the Hamiltonian, the changes 
in the above expression for S in first order of the interac- 
tion are minimal. In addition to the above combination 
g{[p+,p+\ + [p_,p_]) we should add gi[ps,h]^ but since 
[p+,P+\ ^ [P3,A3] ^ [As, As] we obtain 



T.^i{\i + \Q)59/^/2, 



59 



\K{gsui9 + {g + 233) sin 9 cos 9) 



(25) 



This equation (considering it as a precursor to the RG 
equation) is equivalent to Eq.(8) in^^ (see also^). 

In our previous work,^^ we showed that the renormal- 
ization of the impurity in the Luttinger liquid can be 
analyzed within the fermionic formalism. The change of 
9 induced by the interaction and calculated in first order 
perturbation theory in (25) may now be used to obtain 
the renormalization group equation for 6 in lowest order 



[g sin 9 + {g + 2gz) sin 9 cos 9) 



(26) 



In our earlier work we showed that higher order terms 
in the interaction may be summed in a systematic way, 
to access the strong coupling domain. In particular, we 
showed that the one-loop contributions to the RG equa- 
tions for the S-matrix form a ladder series, which can be 
resummed by solving an integral equation of the Wiener- 
Hopf type. The result of this summation reproduces the 
known results obtained with the bosonization method for 
the weak and strong impurity case. It is found to be uni- 
versal in the sense that it does not depend on the choice 
of regularization of the logarithmic divergences in the 
theory. We also showed that two-loop RG contributions 
are absent and the three-loop RG corrections are not uni- 
versal but arc not very sizeable in the whole range for a 
realistic choice of model parameters. 



Let us now discuss how our method can be extended 
to the situation of the Y-junction. First we discuss the 
ladder summation, then we present the results of a com- 
puter symbolic calculation of the three-loop terms. 



2. Ladder summation 

Our previous solution of the ladder equation^® 

amounted to a the dressing of the interaction in the pres- 
ence of the impurity. It thus led to a replacement g ^ g 
in an equation analogous to (26) with 
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1 + ^/rr^-gY 



(27) 



where Y = Tr{p^p^) : notice the different sign in front of 
gY here, which is a result of the different parametrization 
of the S-matrix in our previous work. 

In view of the symmetry 1 O 2, we have Tr (A3. Ag) = 
rr(A3.As) ~ 0. This means that the symmetric den- 
sity combinations remain orthogonal after the scattering 
Tr(p_p+) = Tr{p-p+) = Tr{p-p+) = Tr{p-p+) = 0. 
This property allows us to perform the ladder summation 
separately in each of the channels, as described below. 

First we consider the simpler case of 53 = 0, the ab- 
sence of interaction in the third wire. From Eq. (21) we 
have Tr(p_p_) = cos9 and Tr{p+p+) = \{3 + cos26'). 
The ladder summation of the onc^-loop RG contribu- 
tions in each channel is performed along the previous 
guidelines'^ and the resulting RG equation is written as 



^^^(0) = _£sinfl ^ ^ 
dA ^ ' 4 \l + d- gcost. 



+ 



2 cos 9 



1 + <^- (ff/4)(3-Kcos2e) 



(28) 



with d = 1/ 1 — g2 

Let us now discuss the simultaneous presence oi g, g^ 7^ 
0. The result of the summation in the p- channel is the 
same, whereas the channels p+ and ps begin to mix in 
higher orders of the interaction. We omit the details 
of the derivation below and only provide intuitive argu- 
ments for the result obtained. 

The dressing of the interaction occurs in the individual 
wires, characterized by two density components p+,P3, 
which are orthogonal before the scattering, Tr{p+pz) = 
0. If these density components were orthogonal to those 
after the scattering^ (Tr^p^p^) = Tr{p^p^) = 0, etc.) 
then the dressed interaction would simply be given by 
(27) with y = 0. We define renormalized interaction con- 
stants in the form of two "charges" g, (Z3 by 17 — > 2g/ (1 -|- 

VI-5') = 2?-\ <?3 ^ 233/(1 + ^l-gl) = 2(Z3"' and 
introduce the diagonal matrix Q = diag[q, 93] here. 

According to Eq. (4) the two density components 
/9+, p3 are connected to a vector of scattering eigenmodes. 
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= (Ao,A8)/Aby 



U 



Ao 

±[V2, 
73 I 1, 



(29) 



with U"^ — 1 and Tr{(t)i<pj) = 5ij. Notice that the 
modes Aq and As remain orthogonal both before and 
after the scattering, which is reflected in the matrix 
Y = Tr{(t)i4>j) = rfiag[l, i(l + 3cos26')]. The result 
of the ladder summation is represented by the matrix 
{Q-U.Y.U)-\ 

The rcnormalization of the S'-matrix in the first order 
occurs due to the non-vanishing commutators of densi- 
ties [ps,p3\. Going to higher orders in the in- 
teraction, we have to take into account also the mixed 
commutators, [p+,p3], [p3,p+]. Noticing that the only 
non-commuting components here are Ag and Ag, and us- 
ing Eq. (29), we find that the effect of renormalizing the 
interaction amounts to the replacement 



(5 + 2c/3)sin26' -S^ 6 U.Q.U-Y 



sin 26* 



22 



(30) 



which leads to the explicit result for the /3-function in the 
ladder approximation: 



Q 
q 

K 



2 sin 26* 



+ 



sm( 



^Q — cos 261 q — cos6 
Aqq3 - 2g - 3g3 -H 

2(7 + 93-3 ' 
1 + K I + K3 



(31) 



1-K' 



1-5 



1 



K3 



1 



53 



1 



53 



where we also provide the alternative definition of q,q3 
through the Luttinger parameters K,Kz. 

The equation (31) is a central result of this paper, and 
we analyze it in some detail in the next subsection. It 
should be noted that in the limit 53 ^ and for small g 
this equation reduces to the Eq. (1) of Ref. [22], where the 
fixed points of the point contact model were discussed. 
We confirm Eq. (31) by calculating the perturbative cor- 
rections to the 5-matrix up to the third order of the 
interaction in Sec. V. 



3. Explicit solution of the RG equation 

Let us analyze the RG equation (31). Introducing the 

variable x = cos 0, which determines the conductivity 
components via (22), we rewrite it in the form 



dx 



1/, 2n 42: 1 

2^ ^ \Q + l~2x-^^^^ 

3 2 ^ Ax - X3){x - X4) 

2^"^ 'Ux^-qo){x-q) ' 



(32) 



with 



90 



Xi,i = ^[9 T V3goTg^] 



(33) 



Here Q is defined in (31). The zeros of the r.h.s. , the 
four fixed points (FPs) of the RG equation, are given by 
X\^2 = il) and a;3_4 . For later reference we label the 
fixed points xi,X2,X3,X4 as N,A,M,Q. Fixed point N 
corresponds to Ga = Gb = 0, i.e the three wires are sepa- 
rated. At fixed point A we have Ga = 1, G;, = 0, i.e. the 
third wire is disconnected whereas the main wire is per- 
fectly conducting. The fixed points M, Q arc located on 
the boundary of the physically accessible region, unless 
they are outside the region defined by (20). 

For repulsive interaction, < K^K-^ < 1, the fixed 
point 0:4 > 1 is outside the physical domain, |a;| < 1. 
One can easily verify that the two fixed points N, A at 
xi^2 = ±1 are stable, i.e. dj3/dx\xi 2 < 0- We will see 
below that A turns unstable in the general symmetric 
case. Fixed point M at X3 e (—1,0) is unstable. At 
weak coupling, we write K = 1 — g and K3 = 1 — g^ with 
< 5^ < fif, 53 < 1, and get X3 ~ -g/{g + 2^3), cf.^^. 

The equation (32) may be integrated to give the im- 
plicit solution for x as a function of the length L in the 
form 



F[xiL)] = {l-xy"{l + x) 
= F[x{Lo)]L/Lo, 



■x^r^lx-x^r 



(34) 



where 



71 



72 = 



73 



74 



1 (l-go)(l-g) 
3(l-a;3)(l-ar4) 

1 (l-go)(l + g) 
3(l-hX3)(l-ha;4) 

2 q-X3 qp-xl 

3 a;4 — 0:3 1 — X3' 
2 q- X4 qo- xl 



(35) 



3 a;4 — a;3 1 — x\ 



,2 ' 



We see that the fixed point iV at .t = 1, corresponding 
to three fully detached wires, is characterized by the ex- 
ponent 71, which is determined by the sum of two bound- 
ary exponents, and K^^ . The second fixed point A 
at a; = — 1, corresponding to the ideal wire 1 and 2 and 
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the detached wire 3, relates to the exponent 72, which 
is governed by the boundary exponent of the third wire, 
K^^, and the bulk anomalous dimension of the fermion 
operator {K~^ + K)/2, the latter quantity defining the 
local density of states in tunneling experiments. 

In the important special case of arbitrary but equal 
interaction strength, Ks = K we have Q = (4g — l)/3, 
qo = (2g + l)/3 = X4, and X3 = —1/3. The expressions 
for 73_4 simplify and we have 



73 



74 



1-K 

(2 + if)2' 

K{3 + K)' 



2 i-K 

3 (2+if)2 



{K-^ +K-2), (36) 



For completeness, we also consider the case of attrac- 
tive interaction KjK^ > 1. First of all, the role of the 
fixed points is reversed, so that the fixed points N, A at 
X = ±1 are unstable and the third fixed point M at X3 
becomes stable. Further, we observe that while for repul- 
sive interaction X4 > 1 always lies outside the physical 
region of x, this is different in the attractive case. For 
simplicity, let us first consider the case of equal interac- 
tion in the wires, K = K3. We have X3 = —1/3 and the 
value of X4 = {2q + l)/3 = ^^^^^-^ reaches the physical 
range of x first and coincides with x = —I at K = S. 
At this point the inverse exponents 72^^ = Tj^ = 0. A 
subsequent increase of K takes —1 < X4 < —1/3 inside 
the physical domain. One can verify that the fixed point 
Q at X4 in this latter case is unstable, whereas the fixed 
point A at the edge x = —1 becomes stable. 

Another interesting possibility for attractive interac- 
tion is merging of the points X3 and X4. It happens at 
3go = -q^ in (33), or at 



K. = 2K 



1 + 



{1 + K)'^{2K -1) ' 



and the position of FP in this case is 

X3 = xi = ^il + K)/il- K). 

The last equations show, that in order to have 1x3,4 1 < 1 
one should let K > 2, K3 > 4/3. We return to these 
questions below in Sec. VI, when discussing the "phase 
diagram" of our model. 

B. General Y junction 

Let us now analyze the general case of the ^-matrix, 
Eq. 10. As shown above in Sec. IV A 1, the RG flow to 
first-order in the interaction is determined by two contri- 
butions, which are proportional to [A3, A3] and [AgjAg]. 
We write 



S-\Sr-l= A(a3[A3, A3] + as [As, Ag]) ., 



(37) 



where 0:3 = —g/A^ag = — (g + 2g3)/12 are introduced 
for brevity. Comparing the right-hand side of the above 
expression with the paramctrization (12) of the S-matrix, 
we can find the corresponding change of the parameters 
9, (j), tp in the renormalized scattering matrix Sr- 

After some calculation we find the following set of RG 
equations 



dK 



= (cos tpsinO + sin ip cos 6 sin ^) 



+ 3as cos sin 9 cos 



dK 

dA 



= as 



cos <p sin tp 



— Sag sin (p cos < 



(38) 



sin^ 

3a3 (cos tp sin ^ sin ^ -|- sin tp cos 9) . 



At first sight, there are three equations for three pa- 
rameters. However, one can check that dtp/dA = 
Scoii^ 6d{sm(pta.n9)/dA, irrespective of the interaction, 
so that only two of the equations are independent. 

Indeed, using the equivalent parametrization of the S- 
matrix in the form Eq. (14), we obtain the RG equations 
in the new variables 9,ip,j as : 



d9^ 
dA 
dAp_ 
dA 
dj 
dA 



sin 9{as cos tp + 3as cos 9) 
1 



a3 sin ^ ( 3 cos 
s'mtp 



cos 9 



(39) 



as 



cos 6 



Therefore, we have only two independent RG equations 
(39) for the variables 9 and ip, which define the compo- 
nents of the conductance matrix (19). The third compo- 
nent 7 is not independent and is determined by 9 and tp 
; it does not enter the conductance. 

It is also instructive to express the RG equations di- 
rectly in terms of the components of the conductance 
tensor (18). We introduce 



a = |Tr(A3A3) = cos ^ cos 

6= iTr(A8A8) = i (1 + 3 cos 2^), 



(40) 



so that the conductance matrix is G = diag{^{l 
a), 1(1 - b)). From (39) we find 



da 
dA 
db 
dA 



-2a3(l + b- 20^) - 2as(l - b)a, 
-2{l-b){a3a + as{l + 2b)), 



(41) 



or, directly in terms of conductance 
dGa 



dA 
dGb 
dA 



= as (8C?„(1 - Ga) - iGb) + iasGbil - 2G„), 

= 2asGb{l-2Ga) + QasGb{l-Gb), (42) 
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Comparing Eq. (39) to the above special case, (V' = 
7 = 0), we observe that the arguments, leading to the 
possibility of the ladder summation in Sec. IV A 2, re- 
main valid. Hence we can use the previous result, which 
amounts to the substitution 



1 1 



as 

"8 



2 q — a 
1 1 



1 1 

2q-l + 2Ga' 
1 1 



2Qi-b 
2g + 93 - 3 



2Qi 



1 



(43) 



with q,q3 defined in (31). The equations (41), (42), (43) 
are the main result of this paper. 

Notice that the RG equation for the quantity 



G± — Gb + {2Ga — 1)^ , 

= f(l-6) + a2 = 1 - cos^ 6* sin^ ■0 , 



due to Eq. (42) has a form 



dG± 



«(1-Gi) 



(44) 



(45) 



which, together with (20) and second line in (42), shows 
that the boundaries for the observable conductances for 
free fermions, G± = 1 and Gb = 0, are the RG fixed lines 
in the interacting case. 

Concerning the character of the fixed points in this 
more general situation, there is a qualitative change. In 
the tunneling case, and for repulsive interaction, we had 
two stable FPs , A, N, corresponding to i) one wire de- 
tached {x — —1) and ii) three wires detached (x — 1). 
The unstable FP M (cca) was located between these two 
limiting cases. 

Now we have two independent components of the con- 
ductance. We find only one truly stable FP : N, the case 
of all three wires detached, Ga = Gb — 0. The previous 
stable FP A with one wire detached {Ga = l,Gb = 0) 
transmutes into a saddle point, so it is not truly stable. 
The third (unstable) FP M is at the boundary of the 
region of conductances in the Ga — G;,-plane allowed by 
unitarity and is unstable in both directions. 

The RG flows for repulsive interactions of strength 
g — 0.2,33 = 0.03 are depicted in Fig. 2. It is inter- 
esting to note that the three representative flow trajec- 
tories emanating from fixed point M indicate nonmono- 
tonic behavior of Gb (black) or of Ga (red) as a function 
of, e.g. temperature. As shown in^^ this behavior at 
53 = appears first in second order perturbation the- 
ory in the interaction g. It was demonstrated there, that 
the intermediate M point at 93 = appears due to the 
competition between the bulk zero-bias anomaly and the 
scattering off Friedel oscillations induced by the third 
wire. This underlines the difference in the origin of the 
intermediate point in Ref. [22] and the origin of the in- 
termediate point in this and other studies^'^'^®'^^'^^ with 
93^0- 




FIG. 2: (Color online) Ladder approximation for the /3- 
function and RG flows for g — 0.2 and — 0.03. Three 
fixed points are shown by dots and are : (i) stable one at 
Ga = Gft = 0, (ii) unstable one at Ga ^ 0.8542, Gb ~ 0.4979, 
and (iii) saddle-point type fixed point at Ga — 1, Gb = 0- 
The values of conductances, allowed by unitarity, are shown 
as a shaded region. 



The scaling behavior in the case of attractive interac- 
tion is even more interesting (although attractive effec- 
tive interaction is not easily realized in nature!). The 
conductances at M are given by Gq = (1 — X3)/2 and 



Gft = 1- 



-3' 



where 2:3 is given in (33). As mentioned 



above, fixed point N becomes unstable, while A4 is sta- 
ble. In the case of not too strong attractive interaction, 
1 < K, K3 < 3 , M is the only stable fixed point. Again 
the behavior of both, Ga and Gb may be nonmonotonic. 
For even stronger attractive interaction the fourth fixed 
point Q enters the physically accessible region as a fur- 
ther unstable FP. This leads to a switch of fixed point 
A from unstable to stable. In Fig. 3 the flow diagram is 
shown for strong attractive interaction, K = = 3.7 . 

It is useful to analyze the RG equations (41) around the 
N and A fixed points. For the N point which corresponds 
to a = 6 = 1 we write a=l — a',b=l~b' and find in 
the leading order of the small deviations a' > 0, 6' > 0: 



dA 



{a'-\b')^2{l-K-'){a'-\y), 



dA 



6' ~ (2 - K-^ - K^') b' 



Similarly, around the A with a = 
a = —1 + a' , b = 1 — b' to obtain: 



(46) 
1 we expand 



dA 



(a'- 



i6')^2(l-if) (a'-ifo') 



dA 



b' 



1 - - 



2K 



b'. 



(47) 



In the equations (46), (47) the combination a' — |6' > 
measures the distance to the parabola b = (3a^ — l)/2 
of the simple tunneling case ip = in (40). Thus we 
confirm that the border lines of the physical sector of 
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FIG. 3: (Color online) Ladder approximation for the /3- 
function and RG flows for K — K3 = 3.7. Four fixed points 
are shown by dots and are : (i) unstable one at Ga = Gb = 0, 
(ii) stable one at Ga — I, Gt = 0, (iii) stable one at Ga = 2/3, 
Gi, = 8/9 and (iv) new (appearing a,t K > 3) unstable fixed 
point at Ga — 0.9136, Gb — 0.3159. The values of conduc- 
tances, allowed by unitarity, are shown as a shaded region. 



obtain from Eq. (41) 



da 
dK 
dc 
dK 



(1 - a)(l + 3a) 

q — a 
1 + 3a^ - Adq 



(48) 



{q-df 



This set of equations shows that the RG flows form a 
"spindle" shape around the hnc a = 6, with two edge 
FPs at a = 6 = 1 (N point) and a = 6 = -1/3 (M point). 
The first equation in (48) is easily integrated with the 
result (cf. (34)) 

F[d{L)] = (1 - df- (i + df'' = F[a(io)] L/Lo (49) 

here 7^^ = —2{K^^ — 1) follows from (35) and 

1 - K 



{2 + K) 



(50) 



Notice that the similar-looking 73 in (36) defines the scal- 
ing exponent for c in (48). 



conductances, Eq. (20), are the "fixed lines" of our RG 
equations. 

The meaning of scaling exponents in (46), (47) is as 
follows. The quantity a' — ^b' is the difference of the 
Y-junction from the simple tunneling case which may 
arise primarily due to impurity scattering located in the 
main wire. It is then quite natural that it scales at N 
point with the usual double exponent, 2(1 — K~^), of 
the weak tunneling between two equal Luttinger liquids. 
The quantity b' stands for the tunneling conductance and 
its exponent combines the weak tunneling contributions 
from the main and the third wire. At the fixed point A 
we have the conductance Gb = §6' defined by the weak 
tunneling from the third wire, (1 — K^^), into a perfect 
Luttinger liquid with a bulk exponent (1 — K)'^ /2K. The 
small value of a' — in this case corresponds to a weak 
barrier in the main wire, which scales with the exponent 

2(1 - i^).22.28 



C. Fully symmetric Y junction 

In this section we consider the case of full symmetry 
between the wires, which means equal strength of inter- 
actions, K = K3. The asymmetry at the Y-junction may 
still lead to an asymmetry in the physical properties of 
our system. However, it is a remarkable fact that the 
symmetric Y-junction, a — b, in other terms Gb — ^Ga , 
remains a fixed line in this case. 

We observe that for K — K3 we have q = qs = Qi and 
0^3 = 5(0" Qfs = ^ q)^^ in (43). Introducing 

the symmetrized variables, c — b — a, d={a + b)/2, we 



V. PERTURBATION THEORY AND 
NON-UNIVERSAL TERMS 

The ladder summation discussed above captures con- 
tributions of a certain type (one-loop) to all orders in 
the interaction. By its structure it is similar to the RPA 
summation scheme, which works perfectly well for the 
description of the bulk properties of a Luttinger liquid, 
due to the absence of multi-tail fermionic loops in the 
linearized dispersion model^^. There remains the ques- 
tion of the importance of higher loop contributions. In 
our work on the Luttinger liquid with barrier we were 
able to show in perturbation theory that the two-loop 
contribution vanishes and the three-loop contribution is 
subdominant in the neighborhood of the fixed points. As 
a result, the exponents of the power laws of the conduc- 
tance in length L or temperature T turned out to be 
in exact agreement with those obtained by methods of 
bosonization, and the thermodynamic Bethe ansatz. In 
this section we explore to which extent a similar result is 
true in the case of a Y-junction. 



A. Simple tunneling case 

The symbolic computation of the diagrams up to the 
third order can be performed as explained in the previous 
work.^^ The diagrammatic rules are similar, with the only 
difference that we have 3x3 matrices for the vertices 
now. This means that the kinematic structure for each 
diagram is the same both for the y-junction and for the 
barrier in the Luttinger liquid. Therefore we can use the 
previously obtained results for the individual diagrams, 
while supplying them with different matrix prefactors. 
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We may represent the renormaUzed scattering matrix, 

Sr, in the form 

SKSr = l+9^1+9^^2+9''^3 + ---, (51) 

(in the presence of both g and 173 we consider Ej as de- 
pendent on the ratio 93 /g). In the next step, we verify 
the functional form of the matrix .Sr = expi69{X4 + 
and determine the scalar quantity S6 as a func- 
tion of g, 9, A, 53/5. We have for the renormalizcd value 
Or = 6 -\- 69{9) = f{9,gA,g) and inverting this equality 
we may write 9 = f* {9r,gA, g). Demanding now that the 
initial "bare" value 6 should not depend on A, we may 
derive the RG equation. 

Following these steps, we arrive at an expression for 
the /3-function, which we keep up to third order in 5,53. 
We compare this pcrturbative result with the Taylor ex- 
pansion of Eq. (31) and confirm the identity of the two 
expressions up to second order in the interaction. In the 
third order of g, the direct calculation of the /3-function 
provides the ladder contribution given by the correspond- 
ing term in the Taylor expansion of Eq. (31) and in addi- 
tion a three- loop contribution to j3{9) beyond the ladder 
result. In accordance with Ref. [19], we find that this 
three-loop correction is a sublcading contribution, stem- 
ming from the non-universal parts of the diagrams. As 
was explained in the previous work, performing the eval- 
uation of the diagrams wc find a few generic integrals, 
which contain non-universal numbers 62, 63, 64 which for 
T = are equal to tt^/G, 2 In 2, 2(ln 2-1), respectively.^^ 
The extra term in the /3-function, which should be added 
to (31) is given by 



h{9) = ^ sin^ 9 ib^h + hh + bih) 
/2 = 1 + 8z + 2z2 + + (1 + 253/5)' 
/3 = (1 - zf - Ags/g 
/4 = 2(l + z)(l + z2) 
2;(cos 6) = {1 + 253/5) cos 9 



(52) 



Comparing with the case of the single wire with barrier 
wc sec that the fj are now smooth functions of 9, rather 
than constants. 

Let us consider the influence of ,03 (^) on the fixed point 
values defined by /3l = of the reduced RG equation 
(31), as analyzed in the previous subsection. Near the 
two fixed points N,A, at 9 = 0,7r, corresponding to 
X = ±1, the term (52) is small, proportional to the third 
power of sin 9. This means that neither the positions of 
these two fixed points, nor their scaling exponents 71,2 in 
(35) are affected by P3{9). 

The position of the third fixed point, X3, is defined in 
the ladder approximation by (5 + 253)2:3 +5 = 0, up 
to higher order corrections. The latter condition reads 
z{x3) = — 1 in (52), and we find accordingly that the fj 
at X3 are given by /2 = -4(2 + 53/5)(l - gs/g), fs = 
4(1 — 53/5), /4 = . It follows that the position of the 



third fixed point 2:3 is affected by the /?3 term at 53 7^ 5; 
and the critical exponent 73 is affected in this case, too. 



B. General Y-junction 

We start as described above in Sec. V, but have to 
slightly modify our approach afterwards. It turns out 
that in the case of a general Y-junction the intermediate 
expressions produced in computer calculation show enor- 
mous complexity in the third order. The formulation in 
terms of the 5-matrix followed in the above analysis be- 
comes impractical. 

However, we may find the logarithmic corrections to 
the conductance directly, in the spirit of our treatment 
of the impurity in the Luttinger liquid. Indeed, ulti- 
mately we need to analyze only the partial contributions 
to the two conductances, Eq. (40), rather than the full 
corrections to the S'-matrix, i.e. 3x3 matrix quantities. 
Thus we may keep only these two partial contributions 
from each diagram, which drastically simplifies the cal- 
culation. 

Another simplification occurring in the present anal- 
ysis of the d.c. limit is the absence of vertex correc- 
tions in the diagrams for the conductances.^^ As a re- 
sult, it suflaces to consider only the self-energy parts of 
the Green's fimctions, corresponding to the above Eq. 
(51). The absence of vertex corrections is easily proven 
in case of T = which we consider in this paper. The 
analysis at T 7^ which was undertaken for the impurity 
in the Luttinger liquid, is more involved in case of the 
Y-junction and shall be given elsewhere. 

Below we explain a few points of this analysis in 
more detail. Each diagram contributing, e.g., to the 
component Ga = (1 — a)/2 in (40), stems from the 
generic expression a = |Tr(A3G^^_|_oA3G^^_|_o). Here 
we introduced the retarded (advanced) Green's func- 
tions G^(^) = Ti'd{±t){[t/j{x,t),iP'f{y,0)]) with .x ^ -00, 
y — ^ 00. When Fourier-transformed and taken at zero 
energy, the Green's functions are stripped from the co- 
ordinate dependence and become simply proportional to 
the 5-matrix, G^ = -iS't, = iS, (we set vp = 1). 
The expansion (51) leads to the renormalized quantity 



ao + gai 



5^02 



■g^a3 



(53) 



where ao = ^^'''(AsAs) and 



ai = ^Tr 
a2 = \Tr 
as = \Tr 



^3(^1^3 + ^^3^1) 



A3(S2A3 + SiA3S]]^ + A3S2 



(54) 



A3(S3A3 + S2A3EI + E1A3ST + A3ST) 



and similarly for hr, where the components Sj of the self 

energy have been defined in (51). The advantage of di- 
rectly calculating corrections to , br is most apparent 
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when dealing with the long expressions for S3. The par- 
tial contribution to as has the form 



A3(S3A3 + A3Sj) 

= i?e(Tr(A3A3E3)). 



A3E3A3 



+ h.c. 



The expressions o^, W obtained in this way show both, 

scale-dependent terms in the form of the logarithm of 
A = L/a and in addition scale-independent terms, as is 
explained at length elsewhere. The result of this calcu- 
lation provides expressions for the renormalized values of 
conductance, a,., 5^, in terms of the bare ones, a, b. In the 
next step of the Callan-Symanzik scheme we arc consid- 
ering here, we invert these expansions, and express the 
bare values via the renormalized ones, keeping terms up 
to in the corresponding series. 

Then we require that the bare values be independent 
of the scale in consideration, i.e. the logarithm A, 

da da da dur da dbr , ^ 

and similarly for db/dA. Solving this system of linear 
equations for the quantities dar/dA, dbr/dA, we find: 



d fa, 
dA [br 



f da/dur 
\db/dar 



da/db, 
db/dbr 



d 

dA 



(56) 



This calculation is best done by means of computer alge- 
bra, since the intermediate expressions are quite compli- 
cated. 

We keep the terms of the order of in the final expres- 
sions. A criterion for the correctness of the calculation is 
the absence of any A-dependence on the right-hand side 
of Eq. (56) up to this order. 

When we compare the final expressions, found in this 
direct calculation, with the first terms of the Taylor ex- 
pansion of the corresponding expressions (41), (43), we 
find complete agreement of the universal (regularization 
independent) parts of the /3-functions to third order. In 
addition to these universal contributions of the "one- 
loop" ladder summation (43), we also find non-ladder 
contributions /S^, which first appear in third order and 
explicitly contain the above mentioned regularization- 
dependent coefficients bj . These contributions are rather 
complicated and arc listed in Appendix. 

Comparing these three-loop contributions with the 
simpler case of one impurity in the Luttinger liquid, we 
should make several remarks. First of all, one can ex- 
plicitly check that these contributions satisfy Eq. (45) so 
that the curve G± = 1 remains the fixed line. Second, 
for small g, 53 these terms do not lead to the appearance 
of extra FPs, but the situation with strong interaction 
IflMffal ~ 1 is, strictly speaking, unclear. However, we 
make a plausible conjecture, that all FPs of full RG beta 
function (which contains three loop, four loop etc. con- 
tributions, beyond the re-summed one-loop terms (42), 



(43) ) lie on the borderline of allowed conductances for 
the non-interacting case. 

The next remark concerns the universality of the phase 
diagram, proposed below in Sec. VI on the basis of the 
expression (42), (43). What is the evidence that no new 
FPs appear at strong interaction, and the structure of the 
phase diagram is independent of regularization ? The an- 
swer to this question goes along several lines. From the 
actual form of three-loop contributions we see that regu- 
larization does not change the position of the interaction- 
dependent FP M only in fully symmetric case. However, 
in the latter case the three loop terms are apparently 
unimportant even in the strongly interacting regime, as 
is suggested by comparison of our results with those in^^. 
The bosonization approach by Oshikawa et al. showed 
changes in the character of A fixed point at strong at- 
traction, K = i^3 = 3, and we see the appearance of FP 
Q at this value. The scaling dimensions of leading per- 
turbations, found in Sec. 10.1 and 10.4.2 of^^, agree with 
our one loop formulas (46), (47) for the whole range of 
interaction strength. 

In summary, we again find that the non-universal 
terms in the /3-functions do not influence the behavior 
at the fixed points N,A, which means that the ladder 
summation is sufficient in the case of repulsive interac- 
tions. As for the fixed point M, of importance for at- 
tractive interaction, we find that the non-universal terms 
are relevant and may change the power law exponents in 
principle. We address this question in some detail in the 
next subsection. 



C. Fully symmetric Y junction 

In the important case of equal interactions, 173 = g, 
the position of the third fixed point (M) in Eq. (32) is 
cos 9 = .7;;j = — 1/3 and remains unafii'ected by /33 in (52). 
Expanding around 0:3 we have 



A3 



18 



.9^(x + l/3)(462 - 63 + ^4) + 



(57) 



so that the scaling exponent along the limiting parabola, 
73 in Eq. (36), depends on the non-universal coefficients 
bj, which were found to depend on the regularization of 
the theory. The situation is however more delicate, be- 
cause the scaling exponent 73 in (50) in the direction per- 
pendicular to parabola remains unchanged by the three 
loop contributions. To see that we expand the expres- 
sions for ^3 listed in Appendix B. In terms of variables 
a, c we have the additional three loop contributions to 
Eq. (48) 



da 
dA 
dc 
dA 



3 262 + 64 
^ 16 

a(l — a) 



(l-a)^(l + 3a)^ + 0(5V) 



(262(5 + 7a) + 63(0-1) 



(58) 



+464(1 + 2a)) + 0(5-^ c^) 
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The second equation here at the M fixed point a = — 1/3 
corresponds to (57). The first equation (58) together 
with (48) show that near the fixed points the three loop 
contributions do not change the scaling exponent 73 of 
Eq. (50). 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we presented a theory of charge transport 
through a junction of three quantum wires, modelled by 
Luttinger liquids. We focused on the case of a Y-junction, 
a set-up symmetric with respect to interchanging wires 
1 and 2. We allow for different interaction strengths g 
and 33 in wires 1, 2 and wire 3, respectively. Our method 
employs a purely fermionic representation, which has the 
advantage that the connection to ideal (noninteracting) 
leads is naturally incorporated. The transition from the 
noninteracting leads to the interacting wire is assumed 
to be adiabatic. We find that at zero temperature the 
scattering process is completely described by elastic scat- 
tering (no excitation of real particle-hole pairs). Virtual 
excitations of multi particle-hole pairs are all-important; 
these processes are described in terms of the renormal- 
ized single particle S-matrix. In terms of diagrams for the 
conductances this amounts to the absence of any vertex 
corrections (at T — 0). 

We extended a theory previously applied to a two-wire 
junction to the Y-junction problem. That theory em- 
ploys perturbation theory with respect to the interaction 
in fermionic language (using, however, the concepts of 
current algebra to systematize the bookkeeping) to de- 
rive the renormalization group /3-function for the con- 
ductance. As shown by us in Ref.^^ an RPA type ladder 
summation of an infinite class of terms of perturbation 
theory may be performed to generate all of the known re- 
sults on the scaling behavior of the conductance, power 
law exponents, crossover behavior, and more, for any in- 
teraction strength and any scattering characteristic of the 
barrier. 

In the general time-invariant case the tensor of con- 
ductances features two independent components Ga, Gh- 
These components are confined to an area bounded by a 
curve B in the fundamental domain < G^, G^ < 1. We 
derive the coupled set of RG-cquations for conductances 
in the ladder approximation. It is interesting to note that 
the fixed points of these equations are all located on the 
boundary curve B. In fact a simplified tunneling model 
leads to conductances located on the boundary curve B, 
and allows for an analytical determination of the fixed 
points and the conductances. 

We probe the validity of the ladder summation by eval- 
uating all contributions up to and including third or- 
der (several thousands diagrams). We classify the con- 
tributions into universal (ladder summation) and non- 
universal, with respect to the regularization (finite length 
L or finite temperature T). For repulsive interaction we 
find that the non-universal contributions to the RG-/3- 
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FIG. 4: (Color online) The number and the type of physically 
available fixed points is shown for various strength of inter- 
actions g and g^, as defined by the ladder summation. The 
positions of N and A points are independent of interaction, 
the positions of M and Q points are discussed in text. The 
stable, unstable and saddle-point FPs are denoted by large, 
small and underlined capital letters, respectively. See text for 
additional details. 



functions are subleading in the scaling regime, indicating 
that the ladder summation is fundamentally correct in 
the vicinity of the stable fixed point. For attractive in- 
teraction we find that the non-universal contribution in 
fact changes the location of the fixed points and the val- 
ues of the exponents. 

We find a rich scenario of fixed points. In total there 
are four fixed points A, M, Q, but not all of them are in 
the physically accessible regime. Figure 4 shows the dis- 
tribution of fixed points in the coupling constant g — 53- 
planc. In each regime the stable, unstable and saddle- 
point (unstable) fixed points are indicated by large, small 
and underlined capital letters. For instance, in the case 
of repulsive interactions FP N describes the totally sep- 
arated wires; it is the stable fixed point . FP A becomes 
a stable FP for strong attractive interaction. It stands 
for wire 3 separated from the ideal wire 1 — 2. For any 
attractive interaction with g — g^ the FP M is stable. It 
corresponds to finite conductances in all ways, the value 
depending on the interaction strengths .9,53. We con- 
jecture that M corresponds to the "mystery point" dis- 
cussed in^^ for the totally symmetric junction threaded 
by magnetic flux. The FP Q, finally, is always unstable, 
of the saddle point character, in the limited region where 
it enters the physical domain. 

Most neighboring regions in Fig. 4 differ by one FP 
M, which appears either at the A point or at the A'' 
point, with the corresponding change in the character of 
this latter point. The situation at the interface between 
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nA and nAMq is different, as both "floating" FPs M 
and Q appear at one point on the parabola (see Fig. 3), 
away from A or A''. We recall, that for equal interaction 
strength, g = ga, the Q point appears first at K = K3 = 
3. 

It might be also interesting to note the existence of 
the "tricritical" point between the phases NAM, nA and 
nAMQ, which happens at K = 2, K3 = 4/3, (i.e., at 
g = -3/5, 53 = -7/25) see Sec. IV A 3. This tricritical 
point corresponds to the situation when both M and Q 
points merge with the A point. 

Comparing our findings with previous studies, we ob- 
serve that there is a correspondence between the scaling 
exponents 73 in our Eq. (36) and (l — K^^) for b' in (47), 
and those numerically obtained by fRG method in^^, re- 
spectively 72(3) in Eq. (56) and jf^{3) in Eq. (48) there. 

We confirm that the exponents obtained by our 
method around N point, Eqs. (46), and A point, (47), at 
K3 = K, coincide with those obtained by bosonization 
inl^ Eq. (10.23), and Eqs. (10.106), (10.107), respec- 
tively. At the same time, the exponent {1~K~^) around 
the N point both in our work and Rcf. [18] differs from 
the exponent 1-{3K)-'^ reported in Eq. (2.11) ofRef. [3] 
; we note that in the latter case the RG flow exists even 
in the absence of interaction, K = 1. It was argued in^'^^ 
that the bosonic theory of the Y-j unction has the dual- 
ity property K ^ 3/K, which corresponds to a change 
from the case of totally separated wires (N point) to the 
case of maximally open Y-junction (D point). At the 
latter the conductance exceeds the value allowed by the 
unitarity of the single-particle S-matrix, as argued there 
possibly due to the formation of Cooper pairs at strong 
attraction. The scaling dimension of the leading pertur- 
bation around this hypothetic D point was thus found as 
K/3, in Eq. (10.30) of^^, and as K/9 in Eq. (3.3) of^, the 
latter value evidently arising due to the above additional 
factor 1/3 around the N point. As a result, bosonization 
studies predict a qualitative change in the scaling behav- 
ior of the system at K > 3,^^ or at K > 9.^ Our analysis 
also shows a qualitative change at K ^ > 3, which 
corresponds to the appearance of extra FPs, but not of 
the exotic D type. We stress again, that the RG flows 
in our study always end at the surface of the (generally 
four-dimensional) body describing the conductance ma- 
trix in the absence of interaction. Particularly, the D 
point (which is a = 6 = — 1 or = 1, Gb = 4/3 in our 
notation, see^^) is not a fixed point of Eqs. (41) and the 
RG trajectories do not end there, even if we start from 
outside this body. Our above analysis of the three loop 
RG contributions confirms this picture. We have looked 
for contributions violating the unitarity condition as pro- 
posed in^^, but did not find any. In our formalism such 
contributions would be generated by vertex corrections. 
However, at T = 0, all vertex corrections vanish in the 
d.c. limit. 

In summary, we have derived a renormalization group 

theory description for the two independent conductances 
characterizing a Y-junction of a Luttinger liquid wire 



(1,2) with interaction constant K and a tunneling tip 
Luttinger liquid wire with interaction constant K^. Wc 
summed up infinite classes of contributions in pertur- 
bation theory (ladder approximation) to obtain the RG 
/3-functions. Additional contributions appearing in third 
order were employed to decide whether the result of the 
ladder approximation in the neighborhood of the stable 
fixed points remained unchanged. This was found to be 
the case for repulsive interaction (where arguments can 
be made that all higher order non-ladder terms should 
also be negligible). In the case of attractive interaction 
non-ladder contributions might change the critical behav- 
ior in certain cases. The existence of further fixed points, 
not captured by the ladder approximation, cannot be ex- 
cluded, although it is not very likely. Nonetheless, we 
find it remarkable that our method allows to determine a 
rich scenario of fixed points and RG-flows, including the 
crossover behavior. The corresponding conductances as 
a function of the scaling parameter are readily accessible. 
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Appendix A: Generators of SU(3) group 



For reader's convenience, we list here the matrices Aj 
used in the main part of the paper. The traceless Gell- 
Mann matrices, Aj, with j = 1, ... 8 are the generators 



of the SU(3) group, discussed, e.g., in' 



Ai 



30 





0> 
Ag = ( 1 

1 0; 



As - 



(Al) 




Together with the unit matrix Aq, they have the prop- 
erty Tr[Xj\k] = 2djk, with j,k = 0, . . . 8. After this 
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normalization, the structure of algebra is determined 
by the structure constants fijk according to [A^ , A^] = 
2i^j fjki^i- In the familiar case of SU(2) algebra (which 
is a subalgebra of SU(3), spanned by Ai,2,3), one has 
fjki = Cjki, with totally antisymmetric tensor; it leads to 
simple mnemonic rules. For the present SU(3) case such 
simple rules arc absent, and in most cases we resorted to 
symbolic computer calculations, which are readily done, 
e.g., in Mathematica. 



Appendix B: fUs for the general Y-junction 

The perturbative calculation of the corrections to the 

conductances, as described above, eventually leads to 
RG /3-functions containing one-loop and three-loop con- 
tributions. The one-loop contributions are resummed 
into expressions (43), and the three- loop contributions 
to da/dA, db/dA are found as functions of a, b, Eq. (40). 



^3a=gHF2b2+F3b3+F4b4), 
P3b = 9HF2b2 + F^b3 + Fib^) 



(Bl) 



where 



i=2 = (-432a'' - 56a^(b - 1)k + (-2Ab'^K^ 

864 

+b (39k2 + 405) - 15^2 + 459) - 4aK {2b^K^ 
-Sb^K^ - 96 + + 9) + 3 {4b^K^ - b^ {5k^ 
+27) + b {k^ ~ 81) - 36)) , 

F4 =^ (-108a^ - 18a^(6 - 1)k - 3a^ [b^n^ 
-2b (k^ + 18) + - 36) - aK {2b^K^ 
-36^ (k^ +s) + 9) - 27(6 + 1)^) . 



(B2) 



864 



(6-1) (216a^ + 72a^(6- 1)k; 

+3a (46^ + 6 (k^ - 45) - 5k^ - 27) 
+4(6-l)(26+l)V), 



^3 ^288"^^ " + ^^"^ ~ + ^^""^ ~ 

Fi =^'— ^ (54a-'^ + 9a^(b - 1)k + 3a (2b^K^ 
* 432 ^ V 



-6 (k^ + 9) _ k2 - 9) + (6 - 1)(26 + 1)2^3) , 

(B3) 

and K = 1 + 29,3/9. These expressions are used for the 
analysis of the fixed point M in Sec. VC. In order to 
compare to the result of the simple tunneling case, Eq. 
(52), and to our previous work, it is more convenient to 
go now from the " conductances" a, 6 to the angular quan- 
tities 6,tp and discuss three-loop contributions to dO/dA 
and dip/dA, which we denote by and ^s^, respec- 
tively. 



1^3$ = 9 sin 6 sin i/» cos V' 



262 + 64 
16 



(B4) 



+ ^sin3 0(/262 + /363 + /464), 



/^SV — 9 V' cos 6 



,262-^64 



(B5) 



+ ^ sin Vsin^ e(/^62 + 7^63 + /i64) 

with the functions 

/2 = 4k cos 0{2 cos^ ip + k"^ cos^ 0) 

+ cos tlj{8 cos 2V' - 7 -I- (2 -h cos 26*)) 

/s = cos V'(3 - 2«;(1 + cos 9 cos 7p) + cos^ 9) (B6) 

/4 = 2k cos ^(cos^ V + '^^ cos^ 6 
+ K COS 9 COS V") + 2 cos 3V' 



and 



(B7) 



3cos^ ^-1,9 

/2 = — COS 20-1 

COS0 

-|- 8k; cos V" + 24 cos 9 cos 2tp 

„, 1 — 3cos2 6* /„ „ 9 9 ^ 

= (3 - 2k -h K^ cos^ 9 

■'^ cosfi* ^ 

— 2k cos 9 cos ij)) 

1 _ 3 cos^ 9 

fi = 2 h 12 cos 61 cos 2^- 

cost' 

-I- 2k(1 -I- 3 cos^ 9) cosV' 



The previously considered case of impurity in the Lut- 

tingcr liquid^^ corresponds to setting = in the above 
equations. From the form of (B5) and Eq. (39) one veri- 
fies that the scaling exponents at the RG fixed points N 
{9 = 0, tp = 0) and A (9 = tt, tp = 0) are not modified by 
the presence of /S^e, /S^^, as the latter functions contain 
higher powers of 9, ip in those points. 



Appendix C: Fixed point M in cisymmetric ceise 

In this section we consider the position of the FP M 
in case when all three bulk interaction terms gj in Eq. 
(2) are different. We restrict ourself by the first order in 
Qj, and use the precursor to the RG equation (37) in its 
general form 

3 

S-\Sr-l = -\AY,9j\Pj,Pj]- (CI) 
We parametrize^^ the 5-matrix in T-symmetric case by 

g ^ giA2C/2gjA3(77-V)/2giA5egiA2?/2 _ (-(^2) 
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The conductance matrix in (18) is then given by 

Gaa = I sin^ ^ (1 ~ COS 6 COS tp) + \ cos^ ^ sin^ 9 , 
Gab = Gba = COS ^ Sin'' e, Gbb = sm'e. (C3) 

so that the Eq. (19) is restored at i; = 7r/2. 

Simple calculations show that r.h.s. of (CI) is zero at 
N fixed point, 9 = 0, tp = 0, and at three A points : i) 
6 = IT, tp = 0, ^ = 7r/2, considered above, ii) 9 = w/2, 



^ = 0, iii) 9 = 7r/2, ^ = tt. The position of M point is 
defined by conditions ip = and 

cos9 = -gs^/{g^^+g2^+gs^), 

cose = {92' - 9i^)/{9i^ + 92^ + 253"') , 

at 51 = 52 = we return to the formulas in main part of 
the paper, ^ = 0, ^ = 7r/2, and cos 9 defined by Eq. (25). 
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